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Abstract 

The cocoons surrounding powerful radio sources can be extensive if the jet that feeds 
the cocoon is light and supersonic. They have been shown to remain overpressured 
with respect to the ambient medium for most of the life time of the sources. The 
observed lobes of the radio sources form parts of these extensive cocoons. We show 
that observations of the lobes of giant radio sources allow one to estimate the density 
of the intergalactic medium (IGM) in which the lobes are embedded. We estimate the 
IGM density to be of the order of a few percent of the closure density of the universe. 
We further calculate the radio power of the overpressured cocoon as a function of time 
and the ambient density. 



Subject headings: galaxies: intergalactic medium, jets; cosmology: miscellaneous. 



1 Introduction 



The standard scenario for strong double radio sources involves a 'cocoon' surrounding the 
core and the jet, and consisting of shocked ambient medium and shocked jet material (Scheuer 
1974; Blandford and Rees 1974). It has long been evident that supersonic, low-density jets 
deposit most of their energy in the cocoon, which acts as a "wastebasket" (Scheuer 1974). 
Later studies and numerical simulations concentrated on the properties of the jet and the 
jet/cocoon interface but not much on the evolution of the extensive cocoon. It was also 
the general belief that the cocoons would expand quickly to reach an equilibrium with the 
ambient medium. 

Recently, however, Begelman and Cioffi (1989; hereafter BC) argued that the cocoons 
remain overpressured with respect to the ambient medium for a long time, and that for 



1 



many sources, they have not yet reached a state of equihbrium. They then used this fact to 
address the problem of jet confinement by the ram pressure of the cocoon. Their simple model 
has been numerically verified (Loken et a/. 1992; Cioffi and Blondin 1992). The picture that 
emerges from these studies is that of a overpressured, jet-nourished cocoon, whose length and 
width depend on the balance of the ram pressure of the ambient medium with, respectively, 
the jet's momentum flux and the cocoon pressure. 

We will not be concerned with the problem of confining the jet in the present work, but 
rather the physical conditions in the cocoon and in the ambient medium. Before EC pointed 
out the tendency of the cocoons to remain overpressured, it was generally assumed that 
the pressure in the lobes of giant galaxies (whose lobes lie in the IGM, much beyond the 
hot corona of the host galaxy) also estimated the pressure of the IGM, as the time for the 
cocoons to reach pressure equilibrium with the IGM was thought to be short. That, however, 
does not seem to be the case (see also Subramanyan and Saripalli 1993). Here we ask the 
question whether observations of the lobes can say anything about the physical conditions 
of the IGM. 

We find that following the simple model of BC, one can estimate the density of the 
ambient medium and one can therefore estimate the baryon density of the IGM from the 
observations of the giant radio sources. We also calculate the radio power of the cocoon and 
its dependences on the density of the ambient medium and time. 

We discuss the density of the ambient medium in §3 and the radio power in §4. Through- 
out the paper, we use a Hubble constant oi H — h^obO km s~^ Mpc We also express the 
relevant densities in terms of the closure density of the universe (i.e., = 1 means a density 
of p = {SH"^ / SnG) . We begin with a brief recapitulation of the important aspects of the BC 
model in the next section. 

2 Overpressured Cocoons 

Consider the extensive cocoons produced by light and hypersonic jets. The lightness of the 
jet means that the density of the fluid material in the jet pj is smaller than that of the 
ambient medium Pa, i.e. t] = Pj/pa < 1 (hereafter, the subscripts j,a,c will refer to the jet, 
the ambient medium and the cocoon, respectively). We will not consider heavy or subsonic 
jets. Heavy jets behave like 'bullets' and do not form cocoons. And since the jet luminosity 
is proportional to (M is the Mach number of the head of the jet with respect to the 
sound velocity of the ambient medium), high luminosity of the classical double jet implies 
hypersonic, if not supersonic, jets. 

The jet power Lj is assumed to be constant in time, for simplicity. The velocity of the 
head of the jet Vh is determined by the balance of the thrust of the jet (~ Lj/vj; we assume 
/3j = Vj/c ~ 1, as in BC) spread over the cross sectional area Ah of the bow shock at the 
end of the jet, and the ram pressure of the ambient medium. This yields (eqn. 1 of BC) 




(1) 
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Notice that Ah can be much larger than the cross-section of the jet itself (as in the 'dentist 
drill' scenario (Scheuer 1982)). BC inferred a value of Ah ~ 30 kpc ^ for Cygnus A from 
observations. In their numerical simulation, Ciofii and Blondin (92) found that Ah increases 
with time (roughly, Ah oc t0.4±o.i fiepending on the value of rj), and therefore, Vh decreases 
with time, albeit slowly. In the light of its slow variation in time, and the uncertainties in 
other variables involved in the problem, we will use a constant Ah below. 

The cocoon pressure Pc is obtained simply from the total energy deposited by the jet 
inside the volume of the cocoon Vc = ev(27rr^)//j, where is the half-width of the cocoon at 
the center, Ih = I Vhdt is the length of the jet head and ey is a volume factor depending on 
the shape of the cocoon (see Fig. 1). For a cylindrical cocoon, ey = 1 and for a biconical 
shape, it is ~ 1/3 (Loken et a/.1992). One thus obtains 

where the second equality comes from balancing the cocoon pressure with the ram pressure 
of the ambient medium. 7 is the adiabatic index of the material in the cocoon. This readily 
yields upon integration an expression for the size of the cocoon, as 

^2 ^ {^^y^\^f\^)-'^ . (3) 
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Eqn (1) and (2) can then be combined to give the cocoon pressure, as 

Cioffi and Blondin (1992) argued on the basis of their simulation that initially the cocoon 
energy is not totally thermalized and it is about equally divided between thermal and kinetic 
energy. However, at later times, the fraction of the thermalized energy approaches unity. 

The last equation (eqn (4)) allows us to estimate the time teq for the cocoon to reach 
pressure equilibrium with the ambient medium, after which the cocoon boundary expands 
at the sound speed of the ambient medium. For powerful radio galaxies, with length-scales 
of the order of Mpc, the ambient medium is the IGM. We can write the equilibrium time 
scale as, 

1.34 X 10-(T„,e)-n7^S hli^ L]% f3fAl%{^) "^^^ yr. (5) 

Here, /i ~ 0.6 is the mean molecular weight, = lO^Ta^ K is the temperature of the 
IGM gas, and k is the Boltzmann's constant. We have expressed Lj and Ah in the units of 
gj^g g-i g^j^j 3Q i^pc2 respectively, in accordance with BC. Whereas the earlier studies 
(BC; Cioffi and Blondin 1992, and Loken et a/. 1992) used a value of Ta ~ 10* K, motivated 
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by hot IGM models for the X-ray background, thus obtaining teq ~ 10^ yr, a much lower 
temperature of the IGM seems more reasonable. The constraint from the COBE limit on the 
Compton y parameter {y < 2.5 x 10~^, Mather et a/. 1993) and the anisotropy measurements 
of the microwave background radiation in arcminute scales < 10~^, Subrahmanyan 
et aZ.1993) pose severe problems with hot IGM models for the X-ray background (e.g., Loeb 
1991). Moreover, success of the discrete source models in explaining the X-ray background 
(Zdziarski, Zycki, Krolik 1993) has diminished the motivation for such a hot IGM. Instead, 
models of a photoionized IGM favour temperatures as low as ~ few X 10"^ K (e.g., Miralda- 
Escude and Ostriker 1990). Therefore, eqn (5) suggests a equilibrium time scale that can be 
even larger than the Hubble time (~ 1.3 x lO^^h^Q yr, for a Q = 1 universe). The cocoons 
then stay overpressured with respect to the IGM for all relevant times in the evolution of 
the giant radio sources, unless Lj,Ah have values much different than those used above. 

We now turn our attention to the determination of the ambient density of radio sources, 
in particular, the IGM density from the observations of giant radio sources. 



3 Intergalactic Matter 

The pressure jump across the cocoon boundary can be calculated if the bow shock surround- 
ing the supersonic jet is approximated as a steady, oblique shock (Loken et a/.1992). The gas 

in the cocoon is separated from the shocked ambient gas with pressure Psh,a just inside the 
bow shock (i.e., downstream) by a contact discontinuity (i.e., Pc ~ Psh,a)- The pressure of 
the ambient medium Pa and the cocoon pressure Pc{= Psh,a) are then related as (e.g.. Landau 
and Lifshitz 1959), 

_ 5M^^sin^0-l 
Pa 4 ' 

where Mh is the Mach number of the head of the jet with respect to the sound speed of the 
ambient medium {Cga)- In other words, = v^/csa- 4> is the angle that the bow shock 
makes with respect to the jet axis and 7 = 5/3 is assumed for an ideal gas in the above 
equation. 

For a large pressure jump (pc/Po ^ 1); eqn (6) can be rewritten in terms of the ambient 
density pa, 

3f ^ sm (j) 

Furthermore, since = Ih/t (for a constant Vh), this can also be written in terms of the 
age of the system t and the length of the jet Ih. Note that, eqn (7) is essentially eqn (2) in 
disguise: the transverse growth rate of the cocoon (^) is replaced by a term proportional 
to {vh sin 0) , since the former is hard to estimate. 

(i) Radio galaxies inside clusters: 

Before going over to the case of giant radio galaxies embedded in the IGM, let us use 
eqn (7) to estimate the ambient density for Gygnus A, which is inside a cluster, and whose 
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ambient density has been measured by X-ray observations. This will be a consistency check 
for our method. For Cygnus A, eqn (7) gives the ambient particle density na{— -^) as, 

n„«2xl0-^(g^^^^)(^^^^) (g^) (^^^^) ■ (8) 

From X-ray observations, one estimates the ambient pressure ~ 8 x 10^^^ dyne cm~^ and 
an ambient electron density of ~ 6 x 10~^ cm~^ (Arnaud et a/. 1987). A minimum age of 
~ 3 X lO^Lj^lg yr is estimated from the minimum total energy (~ 10®° ergs) in the radio 
lobes (BC). The pressure jump is, therefore, of the order of ~ 3 if the angle of the bowshock 
is ~ 20°, for which the aspect ratio of the cocoon, i.e. the ratio of the width to the length, 
is ~ 0.4. 

(ii) Giant radio galaxies in the IGM: 

Clearly, the angle is a difficult parameter to measure from the radio observations of the 
lobes. However, for (j) < 30° sin0 ~ tan0 = Vc/lh^ the difference being of the order of unity. 
We can then use eqn (3) to eliminate (j) and rewrite the expression for the ambient density. 
For giant radio galaxies, the ambient density correspond to that of the IGM. We can express 
the IGM density as, 

S1,„M « 0.01A-(1 + {^—^f'\J^y'\,lPfijy^Al%\Zi, . (9) 

We have written the length of the jet Ih in the units of 1 Mpc, Vh in the units of 0.05c where 
c is the velocity of light, and the pressure Pc in cgs units. From the statistics of radio sources, 
Blandford and Rees (1984) estimated the limits of Vh as being 0.03c < Vh < 0.1c. Therefore, 
an average value of 0.05c is reasonable. (Note that Vh is not a constant for all galaxies, as 
evident from eqn (1), and it increases with Lj. We have used an average value in eqn (9) 
only to give a feeling for the numbers. Also, note that there is no real evidence that such a 
value of Vh is also reasonable for giant sources; but here we assume that it is.) The pressure 
in the radio lobes of giant, with Mpc length scales, radio sources at low redshifts {z < 0.1) 
has been estimated from minimum energy arguments, as being of the order of 10~^^ dyne 
cm~^ (Wagget, Warner and Baldwin 1977 obtained pressures ~ 5 x 10^^^ dyne cm~^; but 
see Subrahmanyan and Saripalli 1993). Recently, a giant source 8C 0821+695 ai z = 0.538 
has been discovered with 1^ ~ l-5h^Q Mpc and Pc ^ S x 10^^'^ dyne cm^^ (Lacy et a/. 1993). 

An estimate of the jet luminosity (Lj) can be obtained from the product of the pressure 
in the hotspot, its apparent area and Vj/2 (e.g., Begelman, Blandford and Rees 1974). 

It is worth while to note here that the estimates of the lobe pressure Pc (from minimum 
energy argument) scale as ht^Q^ and Ih scales as h'^Q. The overall dependency of 
estimated from eqn (9) on the Hubble constant is then given by hj^^^'^ . 

Requiring a overpressured cocoon, i.e. demanding p^ > Pa, the lobe pressure observations 
can give an upper limit to the IGM pressure. Using the above estimate of the IGM density, 
one can thus put an upper hmit to the temperature of the IGM gas, Tjgm- For Mpc long 
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radio sources, with lobe pressures ~ 10~^^ dyne cm~^ at the present epoch, an upper hmit 
of TiGM ^ 10^ K is then imphed. 

To summarise this section: we have essentially used the transverse growth rate of the 
cocoons of radio sources to estimate the ambient density (eqns (2, 7, 9)), since the expansion 
of the cocoon in that direction is determined by the balance between the cocoon pressure 
(Pc) and the ram pressure of the ambient medium. We have used the shock jump conditions 
to relate the jet head velocity {vh) and the aspect ratio of the cocoon, to the transverse 
growth rate, leading to eqns (7) and (9). For giant radio sources, this gives an estimate of 
the IGM density as a few percent of the closure density of the universe. 

We will postpone a discussion on the cosmological imphcations of this estimate of the 
IGM density till §5, and calculate the radio luminosity of the overpressured cocoons in the 
next section. We will also discuss in §5 the uncertainties in the estimates of the lobe pressure 
that we have used above. 



4 Radio Luminosity 

The radio power of lobes which are expanding with an equal pressure {pa ~ Pc) or expanding 
freely, with a subsonic velocity, has been calculated by various previous authors (e.g., Eilek 
and Shore 1989). Wc will calculate the radio luminosity in the case of the overpressured 
cocoons discussed above, mainly the way it scales with various parameters. We will follow 
the method outlined in Daly (1994). 

The radio luminosity of the cocoon at frequency u is approximately Lr{i') ~ uPr{u), 
where is the radio power at v. If the rate of synchrotron emission per relativistic electron is 
dE/dt at the frequency v and the number of relativistic electrons with energy E is N{E), then 
the radio luminosity at the frequency u can be approximately written as Lr{v) ~ N{E)dE/dt. 
The energy E = 'yrrieC^, where 7 is the Lorentz factor. For a spectral index a at the radio 
frequencies (i.e., the radio luminosity is proportional to z^^"), N{E) = Ntotij / 'Ji)~'^°' , where 
7; denotes the lower cutoff in energy. Ntot is the total number of relativistic electrons in the 
cocoon, and in the limiting case that the total cocoon energy is in these electrons, can be 
written as Ntot ~ {Ljt/jimeC^). 

The synchrotron emission rate dE/dt ^ 1.6 x 10^^^7^-B^, where the relativistic electron 
spirals around a magnetic field with B± as the component perpendicular to its velocity. 
Assuming equipartition between magnetic and gas energy, the magnetic field in the cocoon 

can be written (from eqn (3)) as B ^ (^^^^^f^) ^ t~^^^. This readily leads to the radio power 
at the frequency 1/, given by 

Prir^) - (4 xl0a)a-c.) (^6")^^^- 'P^'^^' ir-%-t^ . (10) 

This predicts a small increase in the radio power with time for a ~ 1 and with the am- 
bient density pa- Radio lobes in higher density environments would therefore have larger 
total luminosities. One should however remember that there is no direct evidence for the 
equipartition of energy that we have assumed (see below) . 
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Note that Pr above is the total radio power of the cocoon. As emphasized by Loken 
et a/. (1992), the radio luminosity will not be uniform throughout the cocoon. If the last 
reacceleration of the electrons happens at the Mach disk, the backfiow velocity (~ sound 
speed in the cocoon, c^) and the magnetic field B would limit the size of the lobe observed 
at a certain frequency i/. The size of the lobe will be ~ CsTs, where is the synchrotron life 
time of relativistic electrons with energies corresponding to the observed frequency, in the 
presence of magnetic field B. 

However, it is possible that one will observe more extensive cocoons if reacceleration also 
happens elsewhere. It is also possible that if one uses beams that are sensitive enough to 
detect large structures, one could observe limb brightening at the cocoon boundary. This 
would provide support of the BC model for extensive cocoons. 

It is interesting to note here that, at a given redshift, or in other words, at a constant 
Pa, Vh oc Lj^f"^ (eqn (1)), which means that for a constant average v^, the length of the radio 

source Ih Vh oc L^J"^ cx Pr{v)'^ . For a ~ 1, this recovers the correlation Ih oc P^'^ that 
some observers claim to have found (e.g., Oort et a/. 1987). Other studies have also attempted 
to explain this correlation, although with specific models of structure of the ambient matter 
(Gopal-Krishna and Witta 1991). 

5 Discussion 

While using the lobe pressure pc in determining Qjgm as in §3, the assumptions behind the 
estimates of Pc must be clearly borne in mind. Lobe pressure is usually calculated from the 
minimum energy arguments, which gives roughly the cquipartition values. One of course 
does not have any direct evidence that equipartition holds in all cases, although sources 
embedded in clusters seem to yield pressure estimates that tally with X-ray observations. 
Also, the efficiency of energy conversion into relativistic electrons and magnetic field is not 
known. All these contribute to the uncertainty in the lobe pressure estimates. 

It is also important to note that in reality most of the powerful radio galaxies are asym- 
metric and the diagram in Fig. 1 should be treated as being only schematic in nature. 
McCarthy, van Breugel and Kapahi (1991) have examined nonthermal radio emission from 
the plasma and line emissions from the thermal gas in a few of these asymmetric sources 
and suggested that the asymmetry is caused by density inhomogeneity. 

Our estimate of the IGM density (§3) has important cosmological implications. First 
we note that a direct estimation of the density of IGM 0,igm has not been possible in the 
past and only indirect arguments has been put forward to place various limits (e.g., Barcons, 
Fabian and Rees (1991)). Primordial nucleosynthesis model limit the fraction of the baryonic 
density of the universe to fib ~ (0.05 ± 0.01)h^o (Walker et a/.1991). The luminous matter 
in galaxies at z — contribute about flgai ~ 0.007, a similar amount in intracluster gas and 
some galactic or cluster baryonic dark matter may be present ( Barcons, Fabian and Rees 
(1991)). They argued that a value of fljoM ^ 0.01 is therefore reasonable. We note that our 
estimate of flioM from eqn (9) is consistent with these arguments. 
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Clearly, a homogeneous IGM is an ideal case and in reality the universe is probably 
clumpy. The size and magnitude of the dumpiness are still beyond our knowledge, though the 
clumping scale is limited by the upper limits of the anisotropy of the microwave background 
radiation, especially in the arcminute scale (Barcons, Fabian and Rees 1991). It is interesting 
to note that the estimates of lobe pressure in Mpc length radio sources do not vary much. 
Whether this indicates a selection effect or not, needs to be confirmed by looking for giant 
sources with low surface brightnesses (Lacy et a/.1993). In any case, our estimates refer to 
the density in the vicinity of the radio source only and one ought to be cautious in drawing 
conclusions about the global homogeneity of the universe from a handful of giant radio 
sources. 

If the IGM is rather uniform, then, one can constrain the temperature of the IGM gas 
from the limits on Compton y parameter from COBE along with the knowledge of Qjgm (as 
in Barcons, Fabian and Rees 1991, though the limit they used is now outdated by almost 
two orders of magnitude). However, this depends on the temperature and ionization history 
of the IGM, especially at high redshifts. Their argument against a hot IGM model for 
the production of the X-ray background radiation would however be strengthened by our 
estimate of flioM by narrowing the possible range of its value. 

6 Conclusion 

We have outlined a method to estimate the density of the intergalactic medium from ob- 
servations of giant radio sources, in the framework of the model of overpressured cocoons 
surrounding these sources. We estimated an IGM density that is a few percent of the closure 
density of the universe. We also discussed the radio luminosity of the cocoons and found the 
model consistent with the observed correlation I oc P°'^. 
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Figure Caption 

Figure 1 - A schematic diagram of the cocoon surrounding a supersonic jet (see also 
the fig. 1 in BC and Loken et a/. 1993). The jet moves with velocity Vj and the head of 
the jet advances in the ambient medium (here the IGM) with velocity Vh and a with a 
cross sectional area Ah- The bow shock makes an angle 4> with the jet axis and separates the 
ambient medium (with pressure Pa) from the shocked ambient medium (with pressure Psh,a)- 
There is a contact discontinuity between the shocked ambient medium and the cocoon gas 
(with pressure Pc) ■ The width of the cocoon at the centre is Tc- 
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